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ABSTRACT 

We construct general relativistic models of stationary, strongly magnetized neutron 
stars. The magnetic field configuration, obtained by solving the relativistic Grad- 
Shafranov equation, is a generalization of the twisted torus model recently proposed 
in the literature; the stellar deformations induced by the magnetic field are computed 
by solving the perturbed Einstein's equations; stellar matter is modeled using realistic 
equations of state. We find that in these configurations the poloidal field dominates 
over the toroidal field and that, if the magnetic field is sufficiently strong during the 
first phases of the stellar life, it can produce large deformations. 
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1 INTRODUCTION 

After the discovery of the Soft Gamma Repeaters 
and Anomalous X-ray Pulsars (Mazcts et al. 1979 
|Mereghetti fe S tella 1995), a model of these sources 
was proposed according to which they are neutron stars 
with a very strong magnetic field; these magnetars would 
have a surface field as large as ~ 10 14 — 10 15 G, and internal 
fields about ten times larger | |Duncan fc Thompson 1992| ). 
However, a clear picture of the structure, dynamics and 
evolution of magnetars is still missing. For instance we 
do not know whether the toroidal components of the 
field prevail on poloidal ones and how intense they are. 
Consequently, we do not know how large the deformation 
induced by the magnetic field on the star is, an information 
which is essential if one is interested in the gravitational 
wave emission of these sources. A deeper knowledge of 
the structure of strongly magnetized neutron stars would 
also help understanding various astrophysical processes 
involving magnetars (intense activity in the X- and gamma- 
spectra, quasi-periodic oscillations, eventually gamma-ray 
bursts) . 

In a recent paper : ( ":olli et al. 20 09), to be referred to 
as Paper I hereafter, we constructed stationary models of 
non rotating neutron stars endowed with a strong mag- 
netic field, in the framework of General Relativity (GR). 
In these models the poloidal field extends throughout the 
star and in the exterior, whereas the toroidal field is con- 
fined into a torus-shaped region inside the star, where the 
field lines are closed. It is worth reminding that these twisted 
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torus configurations have been found to be a quite gen- 
eral outcome of dynamical simulations of the early evolu- 
tion of magnetized stars, in the framework of Newtonian 
gravity. Furthermore, due to magnetic helicity conserva- 
tion, they appear to be stable on dynamical time-scales 
(|Braithwaite fc Spruit 2004| IBraithwaite fc Nordlund 20061 
|Braithwaite fc S pruit 2006), are not significantly affected by 
rotation ([Yoshida et al. "2 006). and do not depend on the 
initial angle between the rotation and the magnetic axes 
QGeppert fc Rheinhardt 2006 1. 

In Lander & Jones l|2009|) twisted torus configurations 
were studied in Newtonian gravity; the maximal relative 
strenght of the toroidal and poloidal components and the in- 
duced stellar deformation were evaluated using a polytropic 
equation of state (EOS) to model neutron star matter. 

In Paper I, we studied the twisted torus configurations 
in GR using a more realistic EOS. We considered a relation 
between the poloidal and the toroidal components of the 
field which is linear in the flux function, and estimated their 
ratio by determining the configuration of minimal energy 
at fixed magnetic helicity, under the assumption that the 
contribution of the I > 1 multipoles is minimum outside the 
star. 

In the present paper, we reconsider the above assump- 
tions: the higher multipoles contribution is not assumed 
a priori to be minimum outside the star, and we allow 
for a more general parametrization of the relation between 
toroidal and poloidal fields. We determine the configuration 
of minimal energy at fixed magnetic helicity, and evaluate 
the stellar deformations induced by the twisted torus field 
by solving the perturbed Einstein equations including all 
relevant higher order multipoles. 

As in Paper I, the magnetized fluid is described in the 
framework of ideal MHD, which is accurate only in the 



2 R. Ciolfi et al. 



first few hours of the star life, when the crust is still liquid 
and the matter in the core has not yet undergone a phase 
transition to the superfluid state. Since the characteristic 
Alfven time is of the order of ~ 0.01 — 10 s, the magne- 
tized fluid could reach a stationary state while the matter 
is still liquid and not yet supernuicQ The magnetic field in- 
duces quadrupolar deformation on the stellar shape and, as 
we shall later show, for magnetic fields as high as those ob- 
served in magnetars this deformation would be large; when 
the crust forms, it would maintain this deformed shape. 
The magnetic field would subsequently evolve on time-scales 



of the order of 



10' 



10° years due to dissipative ef- 



fects like ohmic decay, ambipolar diffusion and Hall drift 
flGoldreich fc Reisenegger 1992[ |Woods fc Thompson 2006) 
|Pons fc Geppert 2007[ )~ 

The structure of the paper is as follows. In Section [2] 
we generalize the twisted torus magnetic field configurations 
introduced in Paper I by dropping the assumption that the 
contribution from the I > 1 multipoles outside the star is 
minimum, and using a more general parametrization of the 
functional relation between toroidal and poloidal fields. In 
Section [3] we determine the stellar deformations induced by 
the magnetic field. In Section [4] we draw our conclusions. 



2 TWISTED TORUS MAGNETIC FIELD 
CONFIGURATION 

In this Section we briefly describe the formalism and the ba- 
sic equations we solve to determine the twisted torus mag- 
netic field configuration; furthermore, we discuss the modi- 
fications introduced with respect to the analysis carried out 
in Paper I. 



2.1 The model 

We assume that the magnetized star is non-rotating, station- 
ary and axisymmetric. The magnetized fluid is described in 
the framework of ideal MHD, in which the effects of elec- 
trical conductivity are neglected. Furthermore, we assume a 
vacuum exterior. We follow the same notation and conven- 
tions as in Paper I. We treat the magnetic field as a station- 
ary, axisymmetric perturbation of a spherically symmetric 
background with metric 



e u(r) di* + e^ r >dr z + r z {d6 A + sin 2 



Mr) 



(1) 



where v(r) and A(r) are solutions of the unperturbed Ein- 
stein equations; the unperturbed four-velocity is it M = 
(e-/ 2 , 0,0,0). 

To model neutron star matter in the core we use 
the Akmal, Pandharipande, Ravenhall EOS, named APR2 
l|Akmal et al. 199 8). and the Glendenning EOS named 
GNH3 (Glendenning 1985| ; the crust is modeled using a 
standard EOS which accounts for the density-pressure re- 
lation in the crustal region, but not for its elastic properties 
(see Benhar, Ferrari & Gualtieri (|2004p 'l. For a neutron star 



1 We remark that even in presence of a stable stratification of the 
chemical composition, a magnetic field as strong as B > 10 15 is 
still allowed to evolve throughout the star on a dynamical time- 
scale (Thompson & Murray 2001). 



with mass M = 1.4 Mq, the APR2 star has a large com- 
pactness (R — 11.58 km), whereas the compactness of the 
GNH3 star is small (R = 14.19 km). 

As shown in Colaiuda et al. (|2008|) . an appropriate 
gauge choice allows to write the potential A M , in terms of 
which the Maxwell tensor is written (-F M „ = <9„A M — d,j,A v ), 
;'i.s 



A u = 



: E,o,v) ; 



(2) 



the two functions S(r, 9) and ip(r,9) describe the toroidal 
and the poloidal field, respectively. Here is considered as 
a first order quantity, O(B). Furthermore, the ^-component 
of Euler's equation yields —ip, r J r — tp,eJ = 0(B i ). Using 
Maxwell's equations and neglecting higher-order terms one 
finds the integrability condition 

(sin6»S, e ) ie V,r - (sin0E >e ) ir V>,e = , (3) 

which implies that sin#E.e is a function of ip: 

sin0E,0 = 0(y>) = CWi> • (4) 

The function ((ip) = pt(ip)/ip represents the ratio between 
the toroidal and poloidal components of the magnetic field 
and characterizes the kind of field configuration we want to 
model. For instance, configurations with ( = constant have 
been studied in Ioka & Sasaki (|2004p . Colaiuda et al. (2008), 
Haskell et al. ((2008)) . If the space outside the star is assumed 
to be vacuum, the toroidal field, and consequently must 
vanish for r > R, where 7? is the neutron star radius; in this 
case, the choice £ = constant yields an inconsistency, unless 
one assumes that surface currents cancel the toroidal field 
outside the star, or imposes that the constant £ assumes 
very particular values. 

These problems do not arise with the twisted torus 
configurations (see for instance Paper I or Lander & Jones 
(2009), Yoshida et al. (2006)) since the toroidal field is con- 
fined in a region inside the neutron star, and the magnetic 
field is continuous everywhere. For these configurations the 
function P(ip) is continuous, and has the form 

0(VO~ 0(1^1 -1), (5) 
where tp = ip(R,n/2) is the value of function ip, which de- 
scribes the poloidal field, on the stellar surface, and O is 
the Heaviside step function. If j3(ip) satisfies Eq. (0, the 
magnetic field 



0, 



(6) 



for r > R becomes purely poloidal, consistently with the 
assumption of vacuum outside the star. 

To find the field configuration we need to solve the rel- 
ativistic Grad-Shafranov (GS) equation, which follows from 
Euler's and Maxwell's equations (see Paper I for details) and 
has the form 



47T 



4-7T dip 



A' 



1 



4-7J-J.2 



[ip.ee — cot#V,e] 



(p + P)r sin 0[co + citp] 



(7) 



The constants Co, ci characterize the (^-component of the 
current density inside the star, which has the form 

.dp 



a -P^ff + (p + P)r 2 sin 
An dip 



'[c + Cl iP + O(B')]. (8) 
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Figure 1. The function 6M/H m is plotted as a function of a^(R) + a|(i?); on the left £o — 0.61 km on the right the cases £o — 0.65, 
0.61, 0.59, 0.58 and 0.52 km^ 1 are shown together for comparison. 



If we now define ip(r, 9) = sin #a(r, 6),o, expand the function 
a(r, 9) in Legendre polynomials 

oo 

a(r,0) =J2 a i(r)Pi(cos9) (9) 
i=i 

and project the GS equation onto the different harmonic 
components, we find a system of coupled ordinary differen- 
tial equations for the functions a;(r). These equations are 
solved by imposing the following boundary conditions, (i) 
The functions ai(r) have a regular behaviour at the origin; 
an asymptotic expansion of the GS equation shows that this 
implies 

ai{r^ 0) = air l +1 . (10) 

(ii) The functions ai(r) and their derivatives a[(r) are con- 
tinuous across the stellar surface where they match with 
the solutions in vacuum, which are known in an analytical 
form (see Paper I); therefore, the ratios a[/ai computed at 
r = R in terms of the interior and exterior numerical so- 
lutions must coincide, (iii) The overall normalization of the 
field is fixed by requiring the I = 1 component of the mag- 
netic field at the pole to be B po i e = 10 16 G; this corresponds 
to ai(R) = 1.93 x 1CT 2 km. 

Once the form of the function /3(ip) and the number n of 
multipoles we want to include have been assigned, the field 
is determined when we fix n + 2 arbitrary constants: the n 
constants ai and the two constants Co, ci defined in Eq. @. 
n + 1 of them are determined by imposing the boundary con- 
ditions. In Paper I the last constant was fixed assuming that 
the contribution of higher-order multipoles outside the star 
is minimum, i.e. by minimizing the function (5^;>i a ?)/ a ij 
for r > R. In this paper we remove this assumption, thus ex- 
ploring a larger parameter space, and fix the last constant by 
finding the configurations which are energetically favoured 
(this is discussed in detail in the next Section). 

To this purpose, we minimize the total energy at fixed 
magnetic helicity. As shown in Paper I, this is equivalent 
to minimize the ratio 5M/H m , where H m is the magnetic 
helicity 

Hm = ^J d A x^^ s F lS Ap , (11) 
and 8M is the mass-energy increase due to the perturba- 



tion that the magnetic field induces on the spherical star. 
This quantity is determined in terms of the functions which 
characterize the far field limit of the spacetime metric. Fi- 
nally, we compute the ratio of the magnetic energy stored 
in the poloidal field to the total magnetic energy, E p /E m . 
The equations to determine these quantities are given in 
appendix [B] 

The GS system of equations admits two particular 
classes of solutions: the symmetric (with respect to the equa- 
torial plane) solutions, with vanishing even-order compo- 
nents (d2i = 0), and the antisymmetric solutions, with van- 
ishing odd-order components (a^i+i = 0). In Paper I the 
choice of minimizing the / > 1 multipole contribution led 
naturally to symmetric solutions, i.e. those with a-2i = 0. 
It is worth noting that even-order multipoles contribute to 
the energy but not to the magnetic helicity; therefore, any 
solution minimizing energy at fixed magnetic helicity cor- 
responds to a vanishing antisymmetric component. Since in 
this paper we still look for minimal energy configurations, 
we shall consider only symmetric solutions. In addition, as 
in Paper I, we shall restrict to multipoles with / < 5 (the 
relevance of I > 5 multipoles is discussed in Paper I, Section 
5.3). 

2.2 Relative strength of different multipoles 

As explained in the previous section, since we remove the 
condition of minimal contribution from higher order multi- 
poles, the boundary conditions are not sufficient to fix all 
the parameters of the problem and we are left with a free 
arbitrary constant. We choose ci as a "free" parameter. 

In this section we shall choose the function P(if>) as in 
Paper I: 

p{i>) = = <Ro (1^1 - i) ■ e(|^AAI - 1) , (12) 

where £o is a real parameter. In the next section we shall con- 
sider a more general form of /3(ip) compatible with Eq. ©. 
The constant ("o determines the ratio between the ampli- 
tudes of the toroidal and poloidal fields. Thus, we minimize 
the energy with respect to two parameters: ci and £o- We 
proceed as follows. For assigned values of Co: 

• we solve the equations for the a;'s for different values 
of Cl 
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Figure 2. The profiles of the tetrad components of the magnetic field B^(0 = 0), B^(0 = n/2), B^(9 = it/2) are plotted as 
functions of the radial distance normalized to the stellar radius. The left panel refers to the MHM configuration (energy is minimized 
assuming that the contribution of the multipoles higher than I = 1 is minimum for r > R); in this case Co = 0.61 km -1 . The right panel 
refers to the minimal energy configuration ME1, obtained with no assumption on the relative strengths of the different multipoles, and 
for Co = 0.59 km" 1 . 



• we compute SM/H m - the quantity to minimize - for 
the corresponding configurations 

• we compute the quantity y/ a%(R) + a\{R), which repre- 
sents the surface contribution of the multipoles higher than 
I = 1. 

In Fig. [1] we plot the ratio 5M/H m as a function of 
\f a\(R) + al(R). In the left panel we fix Co = 0.61 km -1 . 
In Paper I we showed that, under the assumption of minimal 
contribution of higher order multipoles, i.e. 

<J 'a§ + a\ minimum for r > R , 

(to hereafter, this will be named the Minimum High Multi- 
pole (MHM) condition), the quantity 8M/H m is minimum 
for this value of Co- This MHM configuration corresponds to 
the point A on the curve plotted in Fig. [T] 

Since we now drop the MHM condition, the minimum 
of 5M/Hm occurs for a different value of y/al(R) + a^(R) 
(point B in Fig. [I}, which corresponds to the energetically 
favoured configuration for Co = 0.61 km -1 . For an assigned 
value of Hm, the relative variation of the total energy of the 
configuration B with respect to A is of the order of 13%. 
Fig. [T] refers to a star with EOS APR2. Similar results are 
obtained for the GNH3 star. 

In the right panel of Fig. [T]we plot SM/Hm for selected 
values of Co, and compare the different profiles. We have 
explored the parameter space (fo, \J ' a\{R) + a\(R)) (equiv- 
alent to (Co,ci)), finding that the function SM/H m has a 
m inimum (8M/H m = 0.0817) for Co = 0.59 km -1 and 
\Ja%(R) + a|(_R) = 3.6 x 10~ 3 km. It is worth reminding 
that the I = 1 contribution is ai (R) = 1.93 x 10" 2 km. We 
shall refer to this configuration as the Minimal Energy 1 
(ME1) configuration. 

In Fig. [2] we compare the profiles of the tetrad com- 
ponents of the magnetic field for the MHM and the ME1 
configurations. We see that, whereas for the MHM config- 
uration 5(9) and B( r ) are significantly different from zero 
throughout the star, for the ME1 configuration, obtained 
with no assumption on the relative strengths of the differ- 
ent multipoles for r > R, these field components are strongly 
reduced near the axis. Conversely, the toroidal component 
■B(4>) has a similar behaviour in both configurations. The two 
panels of Fig. [2] illustrate how the magnetic field rearranges 



inside the star when the MHM condition is removed. The 
situation can be explained as follows. The magnetic helicity 
H m can be written as 

H m = -2tt dr {A r ip >e - 4>A r , e )d6 ; (13) 
Jo Jo 

therefore, H m vanishes if either tp = 0, i.e. the poloidal field 
vanishes, or A r — 0, i.e. the toroidal field vanishes. In the 
twisted torus model the toroidal field is zero in the inner 
part of the star, therefore H m receives contributions only 
from the magnetic field in the region where Br^ 7^ 0. Since 
in that region the field components of the MHM and ME1 
configurations are similar, these configurations have nearly 
the same magnetic helicity H m . On the other hand, the en- 
ergy 5M receives contributions from the field components 
throughout the entire star, and these contributions are not 
vanishing in the region where B^\ = 0. When we minimize 
the function 8M/H m in the ME1 configuration, the / > 1 
multipoles, which were kept minimum in the MHM configu- 
ration, do not change H m significantly, but they change 8M, 
and since we require 8M to be minimum, they combine as 
to reduce the field in the inner region of the star. 



2.3 A more general choice of the function (3(ip) 

In this section we construct twisted torus configurations, 
choosing two different forms of the function /?, namely 

/?W = ^Co(!VW'|-ire(|v7V'!-i) , (14) 

(note that a — 1 corresponds to (j!2| l). and 

= -A) - i) CT e(|^| - 1) , (15) 

where /3o is a constant of order O(B). A choice similar to (|14|l 
has been considered by Lander & Jones (|2009|> . who have 
studied the field configurations in a Newtonian framework. 
Although Eqns. (|14[) , (|15[) do not exhaust all possible choices 
of the function /3(ip), they are general enough to capture the 
main features of the stationary twisted torus configurations. 
For P given by Eq. (|14[) the magnetic field components 
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The field configurations are now identified by three param- 
eters: (a, ci,£o) for the choice (|14fl . and (cr, ci,/3o) for the 
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Figure 4. The profiles of the tetrad components of the magnetic 
field are shown for the configuration ME3, corresponding to 
given by Eq. 115H : the values of the parameters are given in Eq. 
CD- 



choice (|15p . As in the previous section, we look for the mini- 
mal energy configuration at fixed magnetic helicity; further- 
more, we compute the ratio of the poloidal magnetic energy 
to the total magnetic energy. We solve the system of GS 
equations for I — 1,3,5 (they are given in appendix |A1 for 
both cases), with the boundary conditions discussed in Sec- 
tion 12.11 For each configuration we compute the magnetic 
helicity H m , the correction to the total energy 8M, and the 
poloidal and toroidal contributions to the magnetic energy 
E m . The equations to determine 8M and E m are given in 
appendix [B] The energetically favoured configurations are 
found by minimizing 8M/H m with respect to the three pa- 
rameters. 

Let us firstly consider the case in which the relation 
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between toroidal and poloidal fields is given by eq. (|14[) . We 
find that the minimal energy configuration (for the APR2 
EOS) corresponds to 

ct = 0.18, Co =0.20 km" 1 , 

\Jal(R) + a? (R) = 3.4 x 10" 3 km . (20) 

We shall refer to this configuration as the ME2 configura- 
tion. In Fig. [3] the configurations with a = 1 (ME1) and 
a = 0.18 (ME2) are compared. In ME2 the magnetic field 
has a slightly different shape: in particular, the toroidal com- 
ponent is larger near the surface of the star, and the exten- 
sion of the toroidal field region along the y axis is smaller. 
The ratio of the poloidal magnetic energy to the total mag- 
netic energy inside the star is 

E p /E m = 0.91 for ME1 
Ep/Em = 0.87 for ME2 . 

If, as in Paper I, we include also the exterior field we find 
Ep/Em = 0.93, 0.90, respectively for the configurations ME1 
and ME2. Furthermore, we find that the minimal energy 
configuration is nearly the configuration with smaller ratio 
Ep/E m , i.e. with larger toroidal component (confirming the 
results of Paper I): thus, the a = 0.18 case also corresponds 
to the minimum value of E p /E m which can be obtained with 
the choice (|14[) . We can conclude that if a is not assumed 
to be 1, we can obtain configurations with a larger toroidal 
contribution, but only by a small amount. 

We now consider the choice (|15[) for the function /3. In 
this case the minimal energy configuration (for the APR2 
EOS) corresponds to 

cr = 0.42, /3 = 9 x 10~ 4 , 

\Ja%{R) + a? a (R) = 3.7 x 10~ 3 km . (21) 

This configuration, which is shown in Fig. 3J will be referred 
to as the ME3 configuration. A comparison with the right 
panel of Fig.|3]shows that the ME2 and ME3 configurations 
are very similar. Inside the star E p /E m — 0.88 and, as in 
the previous case, it is the minimum value which can be 
obtained for this choice of the function /3. 

For the GNH3 EOS, we obtain similar results. The min- 
imal energy configuration is obtained with the choice (|14[) . 
and 

a = 0.30, Co = 0.13 km -1 , 

\Ja%{R) + a 2 5 (R) = 5.1 x 10~ 3 km . (22) 

The ratio of the poloidal energy to the total magnetic energy 
inside the star, for this configuration, is Ep/Em = 0.93. 

We conclude that when we allow for a non-minimal con- 
tribution of the I > 1 multipoles and for a more general 
parametrization of the function /3(ifj), the magnetic field 
changes with respect to the MHM configuration found in 
Paper I as follows: the poloidal field near the axis of the 
star is smaller, and the toroidal field near the stellar surface 
is larger. In all cases the toroidal field never contributes to 
more than ~ 13% of the total magnetic energy stored inside 
the star. 



3 STRUCTURE DEFORMATIONS 

In this section we compute the quadrupole deformation in- 
duced by the magnetic field for the twisted torus configu- 
rations previously obtained. To this purpose we solve Ein- 
stein's equations with a perturbative approach: the magnetic 
field, and the deformations it induces, are considered as per- 
turbations of a static, spherically symmetric background. 
The relevant equations are described in Appendix [C] The 
quantity which is relevant to estimate the gravitational wave 
emission of the deformed star is the quadrupole ellipticity 

e Q = § , (23) 

where Q is the mass-energy quadrupole moment (see eq. 
(|C2|) ). and I the mean value of the star's moment of 
inertia. Indeed, if the star rotates about an axis mis- 
aligned with the symmetry (or magnetic) axis with a wob- 
ble angle a, it emits gravitational waves with amplitude 
(Bonazzola & Gourgoulhon 1996) 

ho =i ^fi 2 /|eQ|sina , (24) 

where Q is the star angular velocity. We remind that we 
normalize the magnetic field by fixing its value at the pole 
as Bp i e = 10 16 G (see Section [2}, and that the quadrupole 
ellipticity scales as Bp ole . 

Furthermore, it is interesting to determine the sign of 
the quadrupole ellipticity because, if eq < (corresponding 
to a prolate shape), a "spin flip" mechanism associated to 
viscous forces may arise, as suggested in Jones (|1975[1 . Cut- 
ler l|2002l) . In this scenario, the angle between the magnetic 
axis and the rotation axis would grow on a dissipation time- 
scale, until they become orthogonal, and this process would 
be associated to a large gravitational wave emission, poten- 
tially detectable by the advanced generation of ground-based 
detectors Virgo and LIGO. 

It is well known that while the poloidal field tends to 
make the star oblate, the toroidal field tends to make it 
prolate. Since in our configurations the poloidal field dom- 
inates over the toroidal one, eq is always positive. There- 
fore, twisted torus configurations are not compatible with 
the spin-flip mechanism, and \sq\ is larger for configura- 
tions in which the toroidal contribution is smaller. We find 
that, for the APR2 equation of state, eq = 3.5 x 10~ 4 and 
eq = 3.7 x 10~ 4 respectively for the ME2 and for the ME3 
configurations. 

It is also interesting to consider twisted torus config- 
urations which do not correspond to minimal energy. In- 
deed, it is not guaranteed that the star sets in the min- 
imal energy configuration before the crust forms and the 
stellar matter become superfluid. Therefore, we have de- 
termined the entire range of possible ellipticities for the 
twisted torus configurations analyzed in this paper; we find 
3.5 x 10~ 4 < eq < 4.8 x 10" 4 for the APR2 EOS, and 
8.1 x 10" 4 < eq < 9.6 x 10" 4 for the GNH3 EOS. As 
discussed above, the larger ellipticities are obtained in the 
purely poloidal limit, whereas the smaller refer to the mini- 
mal energy configurations. We note that, as expected, given 
the mass (1.4 M© in our case), less compact stars (GNH3) 
have larger ellipticities. We also note that the values of eq 
we find for the purely poloidal case are similar to the max- 
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Figure 5. Ellipticities versus /3o for cr = 1, with /3(ip) given by 

imal ellipticity found in Lander & Jones (2009), where a 
polytropic EOS was employed. 

Summarizing, the quadrupole ellipticity eq correspond- 
ing to -Bpoie = 10 16 G would lie in the quite narrow ranges 
(3.5, 4.8) x 10" 4 for the APR2 EOS and (8.1, 9.6) x 10" 4 for 
the GNH3 EOS, i.e. 

^(^tP) 2 xHT 4 > (25) 

with k ~ 4 for the APR2 EOS and k ~ 9 for the GNH3 
EOS. 

As an example of the behaviour of eq when the toroidal 
field contribution changes, in Fig. [5] we plot eq versus the 
parameter /?o for configurations obtained choosing /3(ip) as 
in Eq. (|15[) . assuming that cr — 1 and that the contribution 
of the I > 1 multipoles is fixed by energy minimization. In 
particular, we see that the maximal deformation is given by 
the purely poloidal configuration. 



We remark that the ellipticities given by eq. (|25[) are 
larger than the bounds derived by evaluating the maximal 
strain that the crust can sustain (Ushomi rsky et al. 2000| 
IHaskell et al. 200 6). These bounds do not apply to the case 
we study, because in our case the magnetic field is assumed 
to reach a stationary configuration during the first few sec- 
onds of the neutron star life, when the star is still fluid and 
no crust has formed yet; therefore, if the field is sufficiently 
strong, the deformation it induces can be large, and may 
persist as the star cools down and the crust freezes in a 
non-spherical shape (see also the discussion in Haskell et al. 
(HDDS}, Colaiuda et al. ((2DD5J). 

Recent results from the LIGO- Virgo collaboration 
l|Abbot et al. 20081 lAbbot et al. 2009[) put an upper limit 
on the ellipticity of the Crab pulsar, which should be eq < 
10~ 4 . Our study indicates that strongly magnetized neutron 
stars with ellipticities of this order of magnitude may exist, 
provided this strong deformation was built up before the 
crust was formed, and the magnetic field was sufficiently 
strong. 

In order to further substantiate this scenario, an im- 
portant issue which remains to be clarified is whether the 
twisted torus configurations we find are stable. This issue 
will be the subject of a future investigation. 
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4 CONCLUDING REMARKS 

In this paper we construct relativistic models of non- 
rotating, stationary stars, with a twisted torus magnetic 
field configuration. We extend the work done in Paper I, 
by removing the assumption of minimal contribution from 
multipoles higher than I = 1, by considering more general 
forms of the function P(tf>) which describes the ratio be- 
tween toroidal and poloidal components, and by evaluating 
the deformation that the magnetic field induces on the star. 

We find that the non-minimal contribution of the I > 
1 multipoles, and the more general parametrization of the 
function /3, yield some interesting differences with respect to 
the magnetic field configurations found in Paper I: the new 
configurations have a much smaller poloidal field near the 
symmetry axis, and a larger toroidal field near the stellar 
surface. In any event, the toroidal field never contributes to 
more than ~ 13% of the total magnetic energy stored inside 
the star. 

Since the poloidal field always prevails, in the twisted 
torus configurations the quadrupole ellipticity of the star eq 
is always positive, and its maximum value is obtained in the 
purely poloidal limit. As shown by eq. (|25[) . which summa- 
rizes our results on the stellar deformation, eq depends on 
the equation of state of matter: less compact stars can have 
larger deformations. 



APPENDIX A: GS EQUATIONS 

The harmonic expansion of the GS equations (I17[) . corre- 
sponding to the choice (|14p of the function P(ip), gives (if 
we include the I = 1, 3, 5 components in the expansion) 
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4tt 
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The total energy of the system is E = M + SM, where 
M is the mass of the (spherically symmetric) star with- 
out magnetic field and SM is the contribution induced by 
the magnetic field. It can be determined by considering the 
far field limit of the spacetime metric (|Misncr et al. 19731 
IThorne 1980|l . in terms of the function mo(r) defined in 
equation ()C1[) : 



OT 2 | 4 20 
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a 5 (-315 cos 4 6 + 210 cos 2 6-15) 



sin 2 . (A4) 



The harmonic expansion of the GS equations (|19|l . corre- 
sponding to the choice (115[) of the function P(ip), gives 
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where ip is the same as in (|A4I) . 
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The components of the perturbed Einstein's equations rele- 
vant for the determination of mo{r) give, for the choice (|14|l 
of the function /3(ip), 
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(where <5po is the £ = component of the pressure perturba- 
tion). 

Finally, the magnetic energy is (jStraumann 2 004 
ICiolfi et al.~2009]l 



APPENDIX B: ENERGY AND MAGNETIC 
HELICITY 

The magnetic helicity is given by eq. (|13[) where, if /3(ip) is 
given by eq. (|14p . A r ^ and A T are 

A r ,e = ^V<o(|^|-l) CT 6(|^|-l) , 
sm 6 

A r = e-Cof -^(l^i-l) CT 

xO(|^/V| - l)d6>' . (Bl) 



1 TOO ^ /-7T 

E m = =-I r 2 e~dr sin6B 2 d8. (B4) 

2 Jo Jo 



The above formula can be used to compute the relative 
amount of magnetic energy associated with toroidal and 
poloidal fields. 

The equations for the choice (|15[1 of the function /3 can 
be obtained by the following substitutions: 
V<o -> -0o in (]B1|; V'Co -> /So in ®- 
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APPENDIX C: QUADRUPOLE 
DEFORMATIONS 

The perturbed metric can be written as 
<|Ioka fc Sasaki 20041 IColaiuda et al. 2008P 

ds 2 = -e"(l + 2[ft (r) + fe(r)P 2 (cos 0)])dt 2 

+2[ii(r)Pi(cos6l) + i 2 {r)P 2 (cos 0) + z 3 (r)P 3 (cos 0)J dtdr 

+2 sine UiAp^cos^+^Ap^cos^) 
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The quadrupole ellipticity of the star is defined as eq = 
Q/I EL where Q is the mass-energy quadrupole moment, 
given by the far field limit of the metric (|Clj) 



h 2 (r 



Q/r" 



(C2) 



and I is the mean value of the moment of inertia of the 
star. The value of I can be estimated from the limit Q — > 
of the ratio J/Q in a slowly rotating star model (Q, is the 
angular velocity, J the angular momentum). For M — 1.4 
M we have J = 98.39 km 3 (APR2 EOS) and I = 134.6 
km 3 (GNH3 EOS). 

In order to compute eq we need to solve the following 
system of linearized Einstein equations (here we consider the 
choice (O): 
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2 Note that in Colaiuda et al. 1 12008 I I eq. (81) has a wrong minus 
sign. 
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The integration can be simplified by introducing the auxil- 
iary function 
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This generalizes the variable change adopted in Ioka & 
Sasaki ((2004)) . Colaiuda et al. (|2008|) . With the above sub- 
stitution we are left with two coupled equations 
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where we have used the following relation (arising from 
T% = 0): 



5p 2 = -{p + P)(h2 + jJ [ca + c^ip^e sin 2 o cos Odo 



(C10) 



Eqns. (|C8|I . (|C9|I can be solved using the same procedure 
described in Colaiuda et al. ()2008[) . If we adopt the choice 
(|15fl for the relation between toroidal and poloidal fields, we 
proceed in the same way. In this case the final system of 
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